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Structures, Bourbaki et réalité 


par Louis COUFFIGNAL, 


Administrateur del’ Association Internationale 
de Cybernétique (?) 


Une étude importante de notre collègue M. Mehl a montré que les 
notions de structure, au sens de l’organisation scientifique du 
travail (2) et la notion de structure, au sens des mathématiques 
dites «modernes », sont logiquement équivalentes, les symboles 
constitués par un organigramme, tel que le conçoivent et l’utilisent 
les conseils en organisation, et un graphe, tel que le définit la topo- 
logie, pouvant se transformer l’un dans l’autre par des opérations 
logico-mathématiques. 

C’est là un résultat de grand intérêt pour la cybernétique. I] fait 
entrer, en effet, la cybernétique des relations sociales, telle que la 
concevait Norbert Wiener dans son ouvrage The human use of 
the human beings, dans des théories mathématiques. Comme ces 
théories sont déjà capables de représenter les phénomènes concernés 
par la mécanique, la physique et la biologie (*), la théorie qui, selon 
Wiener, constitue la cybernétique, est tout ou partie des théories 
mathématiques que l’on désigne généralement du nom symbolique 
de Bourbaki (t). 

Ainsi doivent tomber les réserves que l’on faisait encore à la 
conception de Wiener : « We have decided to call the entire field 


(1) Texte de la communication présentée au 3° Congrès International de Cyber- 
nétique, Namur, 11-15 septembre 19061. 

(2) L'organisation scientifique du travail est fréquemment appelée aujourd’hui, 
recherche opérationnelle; théorie de la gestion, théorie de l'administration, etc. 

(5) Voir p. ex. Stanley Jones, La cybernétique de la biologie, Collection Coo 
ion et Cybernétique, Gauthier-Villars, Paris (A paraître) et Henri Laporit, Physio- 
logie humaine, Masson et Cie, Paris, 1959. 

(4) Bourbaki est le symbole de désignation d’un groupe de mathématiciens 
français qui ont pris à tâche de publier un exposé complet des mathématiques 
modernes (Hermann et Cie, Paris). 
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» of control and communication theory, whether in the machine 
» or in the animal, by the name Cybernetics (). » 

Et, ayant moi-même formulé de telles réserves et dit que la théorie 
appelée par Wiener «cybernétique» n’est pas encore faite, je 
reconnais en toute sincérité mon erreur: la théorie que Norbert 
Wiener a appelée cybernétique, c’est du Bourbaki. 

M. Mehl montre également, dans l’article précité, comment un 
graphe est équivalent à un système d’équations logiques. Il en 
résulte que la construction d’un organigramme peut être confiée 
à un mathématicien ou à une machine à calculer. M. Mehl montre 
aussi l’équivalence de ces symbolismes avec celui que j'ai proposé 
moi-même autrefois (?), et qui permet de réaliser effectivement toutes 
les opérations de la logique déductive au moyen d’un simple équipe- 
ment à cartes perforées (et pas seulement de prouver que l'on 
pourrait construire une machine à faire des raisonnements déductifs). 

Ainsi, pour reprendre l’un des cas extrêmes considéré par Norbert 
Wiener, s’il s’agit, dans une entreprise, une administration publique 
ou un État, de faire une réforme de structure, on peut confier le 
travail à n’importe quel mathématicien ayant pratiqué l'algèbre 
bourbakiste ou à n'importe quelle machine a cartes perforées, pas 
nécessairement électronique. Bien entendu, il faudra lui fournir 
toutes les données : les formules d’algébre logique à appliquer, et 
aussi, de façon explicite, les données caractérisant la structure que 
l’on veut obtenir — y compris, par exemple, les services qui seront 
créés pour donner emploi à des fonctionnaires en surnombre, à des 
parlementaires en rupture de réélection ou aux militants de théories 
sociales, métaphysiques, scientifiques ou autres ; il faudra dire qui 
commandera et qui mettra son nom au bas des ordres, si ce seront 
les fonctionnaires qualifiés administrateurs ou ceux qui seront 
qualifiés de gouvernants, dire qui pourra, dans cette structure, 
jouer les Maires du Palais, dire si la compétence sera requise des 
fonctionnaires chargés des contrôles, et bien d’autres détails. 

Sans doute, une telle façon de faire les réformes de structure 
contraindrait-elle à un certain sacrifice moral. Du moins aurait-elle 
l'avantage, quelque élevé que soit le coût de fonctionnement des 


() Norbert WIENER, Cybernetics or control and communication in the animal 
and the machine. Hermann et Cie, Paris, 1948. 

En frangais : « Nous avons décidé d’appeler cybernétique le champ entier de la 
théorie de la commande et des communications dans les machines et les étres vi- 
vants. » 

(?) Louis CouFFIGNAL, La mécanisation de la logique, Mem. Artillerie francaise, 
1952. 
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machines a calculer, de faire une économie sur les avantages de 
titres et de matiére que s’attribuent généralement les réformateurs. 

Mais il ne faudrait pas croire qu’une structure administrative 
ou sociale ainsi fondée sur des théories rigoureuses serait dégagée 
des risques d’échec que courent les structures élaborées par des 
hommes de l’espèce courante, au lieu de machines ou de mathéma- 
ticiens. 

Les théories les plus sévérement rigoureuses, quand on les appli- 
que à la réalité, peuvent engendrer l’absurde ou l’inefficace. 

On en trouvera tout près des exemples typiques. 

La formule 


5e là 


est correctement appliquée dans la langue ordinaire si x est remplacé 
par le nom d’un être particulier et E par la définition d’un ensemble ; 
on dit ainsi, clairement et exactement, 


« Pierre est un homme » 
« Pierre est grand » 


Mais, si l’on dit, avec Carnap, 


« Napoléon est triangulaire » 


on dit une absurdité. 

Autre exemple. M. le Professeur Papy, de l’Université Libre de 
Bruxelles, a exposé récemment, dans une conférence, a Paris, 
comment il avait obtenu de jeunes enfants l’usage correct, spon- 
tané, du symbolisme de l’algèbre moderne. Dans une classe d’enfants 
de 8 ans, on trace au tableau le graphe de la classe : chaque table 
est marquée par une croix et l’on écrit a côté le nom de l'élève qui 
l’occupe. Puis, si une fillette a sa sœur dans la salle, on trace une 
flèche allant de cette fillette jusqu’à sa sœur. Et l’on explore une a 
une toutes les tables, c’est-à-dire toutes les croix. En passant par 
Ernestine, on trace une flèche d’Ernestine à Raymonde, et ainsi 
de suite. En arrivant à Raymonde, on a tracé une flèche de Ray- 
monde vers Ernestine. Celle-ci alors s’est écriée : « Je suis ma sœur. » 
M. Papy, en commentant cette expér'ence, a exprimé son étonne- 
ment et son admiration que cette enfant ait découvert intuitivement 
la réflexivité et la réciprocité de la relation, si subtile, en quoi 
consiste l’état de frère et sœur. Il n’a pas pris garde que cette enfant 
avait dit une absurdité. Et c’est notre tour de nous étonner, sans 


admiration. 
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Cet étonnement s’atténue, toutefois, si l’on remarque que lorsque 
Bourbaki ne trouve pas d’absurdité spontanée dans la question, 
il en fabrique. Pour faire entrer la relation entre droites parallèles 
dans la catégorie des classes d'équivalence, il faut que cette relation 
soit réflexive, c’est-à-dire qu’une droite soit parallèle à elle-même. 
Dès lors, on pose la définition suivante : deux droites sont appelées 
parallèles si elles n’ont aucun point commun ou si elles sont con- 
fondues. Mais désigner par un même qualificatif la propriété, pour 
deux droites, de n’avoir aucun point commun et aussi celle d’avoir 
tous leurs points communs, c’est une absurdité. 

C’est une absurdité de même nature que lancèrent les physiciens 
du siècle dernier lorsqu'ils posérent l'existence de l’éther, en lui 
attribuant, selon les circonstances, des propriétés contradictoires ; 
mais l’existence de cet être de fiction était nécessaire pour faire 
entrer l’optique ondulatoire dans les théories de la mécanique 
rationnelle. 

De tous ces exemples, que l’on pourrait multiplier, on retiendra 
que, sur le plan de l'explication de phénomènes observés, l’absurdité 
nait, en certaines circonstances, de l'application de théories déductives 
à la réalité. 

Une théorie est satisfaisante en tant que telle si elle est exempte 
de contradictions : elle est alors rigoureuse. Cette rigueur, au contact 
du réel, peut devenir de l'absurde. 

Une analyse très simple suffit à démonter le mécanisme de ce 
passage de la rigueur à l’absurdité. 

La première représentation du réel par des symboles abstraits 
consiste à grouper les êtres en classes. Une classe est l’ensemble de 
tous les êtres qui ont un même ensemble de propriétés communes ; 
cet ensemble de propriétés constitue la définition de la classe ; on 
dit encore que cet ensemble de propriétés est caractéristique. On 
connait parfois plusieurs ensembles de propriétés caractéristiques 
d’une même classe ; on dit ces ensembles logiquement équivalents. 
Par exemple, la propriété, pour un quadrilatére, d’avoir ses côtés 
deux à deux parallèles, est la définition de la classe des parallélo- 
grammes ; la propriété que ses diagonales se coupent en leur milieu 
est caractéristique de la même classe ; elle est équivalente à la défi- 
nition. Enfin, la classe est représentée par un symbole qui est géné- 
ralement un substantif. Ce symbole évoque dans l’usage de la langue 
courante, l’image de quelqu'un des éléments de la classe, et, dans 
un raisonnement, la définition de la classe ; c’est en ce sens que 
Pascal a donné comme la clef de la démonstration en mathématique : 
«Remplacer le défini par la définition. » Une classe, pour l’élabo- 
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ration de raisonnements logiques, c’est donc un ensemble de pro- 
priétés. 

Le choix des critéres que l’on retient pour constituer les classes 
d’étres que la nature offre à l’observation constitue une convention 
sociale, c’est-à-dire une convention faite par un groupe d’étres 
humains ; par exemple, les botanistes sont convenus de classer les 
plantes en espéces d’aprés la forme extérieure et la décoration des 
organes de la reproduction ; cette convention n’a pas été faite par 
les ethnologues pour la caractérisation des races humaines. 

Ces notions rappelées, les critères de classification des êtres naturels 
sont des attributs, qualitatifs ou quantitatifs, qu'ont fait apparaître à 
l’homme l'observation et la mesure, ces critères s’appelant dans 
ces derniers cas des grandeurs, tandis que les critères de classifi- 
cation des êtres considérés par le Bourbaki sont les opérations aux- 
quelles peuvent être soumis les éléments d’une même classe. La systé- 
matique bourbakiste n’étant pas fondée sur les mêmes critères que 
la systématique des êtres réels, il n’est pas étonnant que l’on ne 
puisse pas toujours appliquer la première sur la seconde. 

Par exemple, l'opération de réunion de deux ensembles, n’est 
valable, dans les opérations sur les êtres réels où tous les ensembles 
soumis au raisonnement sont des classes, que si ces deux ensembles 
appartiennent à la même classe ; on ne fait, en réunissant dans un 
même parc les poules et les canards, ni un ensemble de poules, ni 
un ensemble de canards ; dans la représentation de la grandeur 
d’un ensemble par un nombre, qui est dès lors un nombre concret, 
on dira que l’on ne peut additionner que des nombres de même espèce. 
La réunion de deux ensembles d'êtres réels n’appartenant pas 
à une même classe devient possible, si on fait abstraction des diffé- 
rences de ces êtres pour ne retenir que leurs caractères communs 
et les situer dans une classe plus générale ; on réunira les poules et 
les canards, ainsi que les oies et les pintades, et même les lapins, 
dans la classe ayant pour nom «la basse-cour ». Bourbaki, ne s’in- 
téressant qu'aux opérations, n’a point scrupule à réunir des ensem- 
bles d’espèces aussi différentes que les droites n’ayant aucun point 
commun et les droites ayant tous leurs points communs, et à en 
faire une même classe, la classe d'équivalence des parallèles. 

Dans cette circonstance, la géométrie, qui par rapport au Bour- 
baki est une science d’étres réels, caractérise cette classe par l’une 
des droites qui la composent prise pour spécimen et par la propriété 
qu'ont les autres droites d’être parallèles à la droite spécimen ; 
la géométrie appelle cette propriété de définition la direction de 


ces droites. 
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Nous avons déja noté que cette convention est nécessaire a 
Bourbaki pour que les étres de cette classe présentent la relation 
de réflexivité. 

Mais la géométrie n’éprouve pas un tel besoin de cette propriété, 
et d’énoncer comme un théorème fondamental qu’une droite possède 
sa propre direction. 

Et nous voici sur la voie de cet autre type d’échec de l’application 
des théories à la réalité, qui est d’être inefficace. 

La réflexivité, en effet, appliquée à des êtres réels, n’est autre 
chose que la tautologie 

«A est A» 
«La lune est la lune » 


La tautologie a fort intrigué les logiciens, depuis l’antiquité. Elle 
n’est pas contradictoire, pas plus que 


« Napoléon est triangulaire » 


et elle n’est pas absurde. 

Sur le plan de l’explication rationnelle des phénomènes, elle est 
parfaitement valable. 

Cependant les logiciens, depuis toujours, en ont écarté l'usage. 
Pourquoi donc ? 

La logique, art de raisonner juste afin de découvrir la vérité, a 
toujours été cultivée pour sa valeur pragmatique, comme moyen 
de faire progresser l’homme dans la connaissance du monde. 

La logique est un moyen d'action. 

De ce point de vue, le raisonnement tautologique est sans valeur, 
parce qu'il est 2neffiicace : c’est pour cela que la logique l’a toujours 
écarté. Il est d’ailleurs écarté de même par un calculateur qui 
transforme une expression et qui tient que, s’il aboutit à la même 
expression, il a perdu son temps ; également par l’organisateur qui, 
en construisant un organigramme, ne manque pas de supprimer les 
circuits bouclés. 

Et l’on voit apparaître la simplicité comme une condition de 
l'efficacité d'un raisonnement (1). 

Dès lors, se révèle par contraste ce caractère du Bourbaki, qu'il 
constitue un langage effroyablement compliqué. Rappelons seule- 
ment cette suite de signes (?), où x et y désignent les longueurs de 


(1 C’est aussi une condition de son élégance. 

(?) Louis CouFFIGNAL, Science, technique et cybernétique, Actes du 2e Congrès 
International de Cybernétique, Namur, 3-10 septembre 1958, Association Inter- 
nationale de Cybernétique, Namur, 1061. 
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segments AB et CD et où chacun des autres signes représente une 
opération particuliére 


Vay = (02)(0u)[Rzu Rux — Czu Ozy (v)(Eux EvzVevn)] 


Cette suite de signes veut dire 


CD'=2AB 


Le Bourbaki, donc, n’est pas fait pour servir. 

Et l’on ne saurait citer un seul problème, parmi les problèmes 
posés par la réalité aux mathématiciens, qui, résolu par le Bourbaki, 
n'ait été déjà résolu auparavant et d’une façon au moins aussi simple. 

Le Bourbaki constitue un jeu de symboles qui s’engendrent les 
uns les autres suivant certaines règles, que les Bourbakistes eux- 
mêmes n'hésitent pas à appeler «les règles du jeu ». Ainsi, le Bour- 
baki est fait pour jouer à Bourbaki. 

Comme certaines de ces régles coincident avec des régles de la 
logique, le Bourbaki peut aussi étre considéré comme un forma- 
lisme représentant les opérations de la logique. Mais on a connu 
déja une logique formelle, qui, au Moyen Age, était imposée a tous 
les étudiants par les Docteurs en Sorbonne, et qui a disparu par 
suite de son inutilité. 

Pour autant que l’histoire déroule dans le temps des événements 
analogues, c’est une nouvelle scholastique qui se crée sous nos yeux 
pour exiger des étudiants qu’ils apprennent, généralement contre 
leur gré et sans le comprendre, l’inutile Bourbaki, comme, cinq 
siècles en arrière, on apprenait par force l’inutile Baralipton. 


L’inefficacité des théories qui constituent l'algèbre moderne, 
le Bourbaki, est observable dans presque tous les cas où l’on veut 
l'appliquer au réel. Pour des théories exprimées dans un autre 
langage, les échecs sont plus rares, mais ils existent et demandent 

’être prudent. 

D'abord, le philosophe anglais Thomas Greenwood (1) a montré 
que les opérations de la logique moderne et les opérations de la 
logique aristotélicienne sont les mêmes ; que l’on formalise ses 


() Thomas GREENWOOD, Les fondements de la logique symbolique. Hermann et Cie, 
Paris 1037: 
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raisonnements dans le langage ordinaire ou dans un langage arti- 
ficiel quelconque (#), on risque donc, les mêmes ennuis. 

Mais raisonner, dans la vie courante et dans l’étude des sciences 
selon les méthodes classiques, s’apprend en faisant des raisonne- 
ments ; lorsqu'un raisonnement est absurde et inefficace, on apprend 
par l’action pédagogique (?) de l’environnement que ce raisonne- 
ment ne doit pas être répété. Il s'établit ainsi dans l’esprit de chacun, 
une connaissance plus ou moins consciente du domaine d'efficacité 
des raisonnements que l’on a déjà rencontrés. L’application de ces 
raisonnements, c’est-à-dire de théories, à des cas nouveaux se fait 
par analogie avec des cas déja résolus : les échecs deviennent ainsi 
beaucoup plus rares. Lorsqu'il s’agit de l’organisation de l’action, 
ils tiennent généralement à des causes faciles à caractériser : 

— s'il s’agit d'action sur la matière inerte, l’échec est dû le plus 
souvent à la méconnaissance de conditions technologiques, 

— s’il s’agit d'action sur les êtres vivants, l'échec est dû a une 
réaction imprévue du milieu sur lequel on exerce cette action, 
parce que la réaction d’un tel milieu est du type aléatoire (3), ce 
qu’on exprime encore en disant qu’il jouit d’une certaine liberté. 

Le premier cas s'illustre abondamment dans l’art de l'ingénieur, 
où l’on voit nombre de projets ne point aboutir parce que techni- 
quement irréalisables, et des matériels construits par des ingénieurs 
insuffisamment instruits de la technologie subir des incidents de 
fonctionnements nombreux ou provoquer des accidents et parfois 
des catastrophes. 

Dans le second cas, on trouve ces innombrables exemples d’un 
règlement inappliqué parce que inadéquat aux circonstances ou 
des règlements appliqués de façon rigide qui créent des situations 
ridicules — le ridicule est, dans l’action sur l’homme, un manque 
d'efficacité —- tenant généralement à une sorte de conviction du 
législateur ou de l'organisateur qu’un organisme fait d’étres hu- 


(1) Les théories de la chimie organique, par exemple, sont représentées, selon 
des conventions admises par tous les chimistes, dans un langage qui ne doit rien 
aux mathématiques. 

(2) Le terme est pris au sens large, par exemple celui qu’il a dans la langue 
allemande, où la « pédagogie » désigne les méthodes et les techniques de transmission 
d'informations à des hommes en vue de fixer ces informations dans leur mémoire. 

(*) Selon la terminologie que nous avons proposée: Louis COUFFIGNAL, Les 
notions de base, Collection Information et Cybernétique, Gauthier-Villars, Paris, 
1958 et Essai d’une définition générale de la cybernétique, Actes du rer Congrès 
International de Cybernétique, Namur, 26-29 juin 1956, Association Internationale 
de Cybernétique, Namur, 1958. 
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mains est une machine, alors qu'il n’est qu'un mécanisme, et qu’il 
obéira desmodromiquement aux ordres qui lui sont donnés. 

Dès lors, pour l’homme d’action, une théorie, ou un système de 
principes ou d’axiomes, n’est qu’un modèle qui suggère une solution 
possible au problème qu’on lui fait simuler. 

L'efficacité de ce modèle sera d’autant plus grande, comme pour 
tout autre modèle, 
que la représentation des données par le modèle sera plus 
fidèle et plus complète, 

— que l'extension au cas nouveau étudié du domaine d'efficacité 
connu de ce modèle sera moins lointaine, 

— que le modèle sera plus simple. 

Sous ces remarques générales, l’homme qui a acquis par la pra- 
tique une bonne connaissance du domaine d'efficacité des théories 
n’aura point de peine à choisir celles qui l’aideront le mieux a 
résoudre les cas particuliers qui se présentent à lui. 


CYBERNETICA (Namur) 
Vol. IV — N° 4 — 1061 


Computers and Simulation 


by Carl HAMMER, 


Manager, Scientific Applications, Electronic Data Processing 
Division, Radio Corporation of America (Washington) 


INTRODUCTION 


Electronic computers date back to the early 40’s and we may 
safely assume that by now most everyone is familiar with their 
theory and operations. This presentation will discuss the applica- 
tion of electronic computers to simulation, in particular to the si- 
mulation of communications networks. 

To set the stage, we shall briefly discuss computer theory, some 
hardware aspects, and those coding and programming techniques 
which may appear relevant to the design and development of 
simulation programs on large-scale electronic digital computers. 

During the course of this discussion, we shall also give an example 
for the simulation of a simulation ; not having a computer avai- 
lable, we have desk-simulated the very functions that would have 
been executed on an electronic computer. The subject of our simu- 
lation is, not too surprisingly, a somewhat simplified communica- 
tions network. Its results will hopefully give a better insight into 
the operating characteristics of such networks and into simulation 
techniques in general. 


THE ROLE OF COMPUTERS IN TODAY s WORLD 


Broadly speaking, computers are in use today because they can 
cope efficiently with four functional aspects : (i) problem solving, 
(ii) data processing, (ili) control, (iv) simulation. 

One of the oldest applications of electronic computers is to the 
area of (mathematical) problem solving. We include here the tabu- 
lation of mathematical functions, not otherwise obtainable, or 
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of special data, of interest only to a limited clientele. As is well 
known, the first large-scale electronic computer, ENIAC, was 
designed and used for the calculation of ballistic firing tables. Its 
eminently qualified staff learned quickly that it could also be used 
to calculate 7 to ‘* only ’’ 2000 decimals [13]. This feat of computer 
acrobatics should neither be a surprise nor an indication of the 
computer’s limitations. 

Historically, computers were then successfully employed to 
reduce man’s data processing workload. One of the earliest large- 
scale general-purpose computers to go into operation was the 
UNIVAC, at the Bureau of the Census. Grinding out statistical 
tabulations at an hitherto unknown rate, this installation caused 
also the first, vexing engineering problem which led to the deve- 
lopment of high-speed printers. Visualize, if you will, the frustra- 
tions of the statisticians when, after having processed the data, 
they were unable to print them out ! The advent of the HSP was a 
significant step in the reduction of interface problems between 
computer and man. However, since that time other problems have 
arisen and today’s mountains of paper forms, punched cards, and 
magnetic tapes testify to the need for another ‘* breakthrough ”’. 

More recently, computers have been introduced to the task of 
simulation. Two significant characteristics of simulation problems 
are the mathematical model and its validation. Simulation, as we 
shall see later, is the art of make-believe. The results obtained 
from such programs cannot be better than the mathematical or 
statistical models upon which they are predicated. 

Finally, computers have also been assigned the functions of 
real-time control. As such, they may be destined to assume a signi- 
ficant part of the decision-making functions hitherto reserved for 
the judgment of man. The speed with which such decisions must 
now be made exceeds man’s capacity for reaction by factors of 
thousands, or even millions. Thus computers will speed up decision- 
making processes to effect tighter and more absolute control over 
sensitive parts of large engineering systems. For example, literally 
dozens of computers are at work during the firing of one single 
missile shot ! Their main function is that of control over the flight 
path, from launch through impact. 


COMPUTERS AND THEIR LANGUAGE 


Let us now briefly review what makes computers tick. Specifi- 
cally, let us talk about three items which must be considered above 
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all: hardware, software, and support. Most everyone thinks of 
computers as hardware. A computer installation is viewed as a 
large air-conditioned area, filled with towering racks of electronic 
equipment, which are only overshadowed by mountains of tapes 
and paper. To many users, the hardware is still the focal point of 
interest. In order to select a computer, they will look at the central 
processor, its speed, and at its logic. They will make attempts 
to compare its functions with those of other computers. The same 
type of user may also be interested in the available storage devices, 
core memories, drums, and tapes. Further, he may look at input 
and output devices, such as card readers and punches ; paper and 
magnetic tape units ; photo-electric scanning devices ; and finally, 
high-speed printers and cathode-ray-tube displays. 

The last years, however, have brought about a change. The 
computer industry has begun to recognize that the user is not 
served by hardware alone. A new term has come into being : soft- 
ware. This term describes the library of programs, the ‘* bag of 
tricks ’’, which make the computer perform. Behind this software 
stand the untold hours of toil by programmers, coders, and ope- 
rators ; by the human beings whose efforts will benefit other human 
beings. Special recognition should be given here to the ‘ brain 
behind the brain ’’, the system analyst. It is through his support 
efforts that we make efficient use of these machines, which, after 
all, cost hundreds of dollars per hour to operate ! 

In summary, we see that the application of electronic computers 
is most certainly not simply a matter of setting up two million 
dollars worth of hardware in the hope that nature will take care 
of later events. On the contrary, the successful installation requires 
planning, an abundance of analytical skills and, above all, a devo- 
tion which is not diminished by long and late hours. 


THE CONCEPT OF SIMULATION 


Simulation, as was mentioned earlier, is a matter of make-believe. 
Here is what Webster has to say : ‘ Simulate : to assume or have 
the mere appearance of, without the reality ; to assume the signs 
or indications of, falsely ; to counterfeit ; feign, imitate. ”’ 

The art of simulation is probably as old as little children making 
mud pies. Today the tools used are more sophisticated and costly ; 
hopefully, the results are also more digestible. The primary tool 
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employed in the art of simulation is the mathematical model. 
Because of its importance, we shall give it more than a cursory look. 
We shall also discuss such concepts as system analysis, formulation 
of mathematical models, their evaluation, and machine methods 
used to implement them. 

First, let us propose a general thesis, namely, that everything 
done with computers is merely simulation. Computers do not really 
provide transcendental or irrational functions ; indeed, they provide 
only approximations thereof. Notable along this trend of thought 
are the efforts of Cecil Hastings on the beaches of Waikiki [5], or 
of George Reitwiesner [13] who computed the first two-thousand 
decimal places of 7. But even this enormous stream of digits comes 
not much closer to the ‘* truth ”’ than 22/7, considering 7’s infinity 
of digits. 


MATHEMATICAL MODELS (1) 


The word ‘* model ”’ is generally used in scientific work to denote 
a miniature replica used for the demonstration of ideas or methods. 
In a broader sense, the word may be used to denote the forms, 
patterns, or equations appropriate for the reproduction of indi- 
vidual ideas and for the interrelationships existing within a group 
of ideas. The making of models is among the oldest pursuits of 
mankind. The simpler ones are direct reproductions of obvious 
features of physical objects. The more abstract and sophisticated 
ones bear scarcely any resemblance to the entities represented. In 
fact, models are now often used to represent not only physical 
objects, but also the most abstract notions. 

The systems engineer is concerned with the rigorous manipu- 
lation of individual ideas comprising a group. In order to clarify 
existing or suspected interrelations between them, he tends to 
consider all models as essentially mathematical in character. 

Mathematical models ordinarily can be treated in accordance 
with expressed rules and prescribed logic. Thus, they engender 
confidence in the conclusions without reference to their ‘+ truth ”’. 
Other than mathematical models may be subject to emotional 
interpretation and psychological mood. 

One category of mathematical models consists of the so-called 


# (4) This section and the five following ones were adapted from some unpublished 
material by C. W. Steeg, RCA, Waltham, Mass. His permission to use this material 
is gratefully acknowledged. 
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direct analogs. Most direct analogs are scaled replicas of the mo- 
deled system. The concept of scale models, together with other 
equivalent or analog forms of reproduction, has facilitated the 
performance of many engineering tasks. Often it is virtually im- 
possible to make direct studies on the system under consideration. 
Yet much valuable information can be gleaned from the study of a 
purely scale model of the system. The use of scale models in the 
study of dams in controlling the runoff in a water-fed area is a 
typical example. Wind tunnels for the collection of aerodynamic 
data on proposed air frames represent somewhat more advanced 
applications of the ideas of scale models. 

Network analyzers are another form of direct analog, or scale 
model. In this type of model an electrical power system consisting 
of generators, transmission lines, and electrical loads is represented 
by a collection of lumped-resistance capacities and inductances 
through lumped parameters and distributed portions. Voltages 
representing the generators are impressed on the analyzer ; currents 
and voltages are measured at other points in the network as they 
would be measured in the system under study. 

The demands placed upon modern control systems have taxed 
severely the ability of designers to analyze existing systems or to 
synthesize new ones. Troubles arise especially where one group 
designs one portion of a system, another group designs some other 
portion, and still a third group attempts to unite these two portions 
into an integrated whole. In the past few years the realization has 
been growing that a complex system must be considered as a whole 
if it is to perform effectively. Consequently, a new branch of ‘* sys- 
tems engineering ”’ is emerging. Computers are playing a vital rôle 
in this work. The mathematical models formulated by systems 
engineers are effective only if they are developed correctly and if 
they lend themselves to study by the techniques that the analysts 
can bring to bear for this purpose. The primary goal that the analyst 
should maintain while formulating a mathematical model for a 
system is to develop a model that is sufficiently complete, consis- 
tent, and accurate so that its study will enable him to answer 
questions posed about system performance. As a secondary objec- 
tive, he should strive to develop a model that is not unnecessarily 
complicated because the simpler the model, the easier the subse- 
quent task of analyzing it. Consequently, the analyst should be 
familiar with the techniques that have found extensive application 
as aids either in the simplification of models or in the generalization 
of the results obtained from the study of models. 
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SYSTEM ANALYSIS 


One of the new and vital tools of the systems engineer is simu- 
lation. This concept involves the manipulation of a dynamic model 
and its examination on a large-scale computer, either analog or 
digital. When defined in this way, simulation denotes the study of 
an abstraction of a system oriented toward a specific purpose and 
conveys no relation to duplications, miniature copies, or static 
mock-ups of sets of equipments. Simulation is, however, only one 
aspect of the rôle that computing equipments are playing in systems 
engineering. Therefore, we shall digress briefly to discuss the more 
general concepts of control systems engineering. 

Because the systems engineer is confronted with an extremely 
wide range of problems, no definite set of steps can be presented 
that lead to the solution of any specific problem proposed. However, 
experience has taught us the principal steps through which the 
study of all complex systems must go. 

The most important stages in the analysis of a complex system 
are the following : 


(i) Preliminary study, 
(ii) Formulation of a mathematical model of the system, 
iii) Study of the model, 
v) Evaluation of the results, 
(v) Re-examination of the entire problem, and if necessary repe- 
tition of the four preceding steps: feedback, 
(vi) Construction and test of a prototype system. 


(i 
(1 


Although the boundaries between these stages are often some- 
what fuzzy, this division will prove helpful. 

Before an effective detailed analysis of any new system can be 
made, a broad preliminary study should cover the entire system. 
This preliminary study should lead to a very clear understanding of 


(i) The basic functions that the system is to perform, 
(ii) The basic physical phenomena involved in the operation of 
the system, 
(iii) The form of the system that will be acceptable to the ultimate 
user and that will be compatible with the conditions imposed by 
such factors as environment and economies. 


In a preliminary study every effort must be made to differentiate 
clearly between those characteristics that are essential and those 
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that are incidental. Unless the engineer engaged in the study of a 
complex system is quickly able to pass over the trivial aspects of 
the problem and concentrate on those of fundamental significance, 
many pitfalls await him. 


FORMULATION OF MATHEMATICAL MODELS 


During the next stage in the study of a complex system, the 
analyst should attempt to set up a mathematical model of the 
system. In the study of a complex system the analyst has great 
freedom in the choice of mathematical models, and he may employ 
several different models during the course of the study (See the 
section *‘ Simulation methodology ’’). 

For many systems with which the analyst is concerned, present- 
day knowledge allows description of the phenomena involved in 
such detail that the model for even a very simple system becomes 
extremely complex. Consequently, if a useful model is to be obtai- 
ned, it must be found by ignoring those properties and conditions 
of the system that are irrelevant to the problem at hand. The model 
must describe accurately, or with suitable approximations, only 
those parts of the system that are essential. A model satisfactory 
for a particular investigation depends on the fundamental purpose 
of that investigation. The model chosen may be considered satis- 
factory if good agreement is obtained between the results gained 
from the study of the model and from appropriate tests on the 
actual physical system. During the study of a complex system, a 
variety of models with varying degrees of complexity may be 
employed. In some instances a logical development toward a more 
and more detailed model can be made ; but in other cases, it may 
prove necessary to abandon a particular system, revise some of the 
basic concepts employed, and on the basis of these new concepts 
develop an entirely new system. 

A mathematical model that includes the detail necessary to 
represent a complete system may become very unwieldy. Conse- 
quently, it does not lend itself readily to solution by either analysis 
or machine methods. A study utilizing a complex model is inadvi- 
sable until simpler models have been examined extensively and a 
system appearing basically satisfactory has been obtained. Even 
at this stage tests on the complex model may indicate the need for 
the revision of the basic system. The chance that such will occur 
is decreased enormously if a thorough study of the simpler models 
has been completed. 
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STUDY OF MATHEMATICAL MODELS, 


The difficulties in specifying detailed methods appropriate for 
the solution of the equations used to describe a complete system are 
as great as those in giving definite rules for formulating the equa- 
tions. However, if during the course of the study the complexity 
of the models employed has progressed from fairly simple cases to 
a more complete representation, the analyst, while developing the 
models, also has had an opportunity to develop methods of solution 
appropriate to the particular problem. 

If practical analytic solutions can be obtained, the analyst is in a 
much better position to predict what effect the various parameters 
in the system have upon system performance and to arrive at a 
system that is optimum in some sense. However, the fact that an 
analytic solution to a problem can be obtained does not necessarily 
mean that the system engineer will consider his problem solved and 
have no need for machine computation. The analytic solution to a 
problem does not necessarily take the form most desirable for 
design purposes. The mathematician may feel his task is completed 
even if the solution is given in terms of some integrals that cannot 
be evaluated analytically. Such a solution is often of little value 
to a design engineer who would consider much more valuable the 
results of numerical calculations for a few well-chosen cases. 


EVALUATION OF MATHEMATICAL MODELS 


In certain types of problems the task of the analyst is essentially 
completed at the end of the machine-computation phase of the 
study. However, in problems presenting a real technical challenge, 
the analyst frequently makes his most significant contributions 
while evaluating the machine results, by reducing these results to 
such a form that they are useful to those responsible for making 
broad design decisions. 

Evaluation of the results of a study is accomplished in several 
steps, which for convenience may be called : (i) analysis of data, 
(ii) formulation of conclusions, (iii) presentation of results. 

An important part of the analysis of the data obtained in a study 
customarily takes place during the time a problem is set up on a 
computer. However, much work usually remains even after the 
computer phase of the study is completed. Only then can the ana- 
lyst place all the solutions obtained in proper perspective and 
reduce the significant results to the form of curves and tables which 
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will be helpful to him in formulating conclusions. Once the data have 
been reduced, the analyst may be able to reach the appropriate 
conclusions immediately. However, frequently it becomes apparent 
that additional work must be done before one can obtain or present 
conclusions of the study. It may prove feasible to obtain additional 
data either from the study of a simplified system or by purely 
analytic means. Possibly, based on the experience gained thus far 
in the study, the analyst can justify the omission of some of the 
simulated effects. Their importance may not have been apparent 
before the machine solutions were obtained and he may thus realize 
appreciable simplifications in the analysis. If simplifications of this 
type can be made, they will reduce the difficulties in obtaining 
additional machine solutions and increase the probability of analy- 
tic solutions. Alternatively, the analyst may have observed on the 
computer the forms of the solutions obtained and the effect of the 
important parameters on the performance of the system. Thus he 
is sometimes able to obtain analytic or quasi-analytic solutions. 
If neither of these possibilities is practicable, further machine study 
with the same setup may prove necessary. 

With simple systems, formulation of the conclusions of a study is 
relatively easy. But with complex systems the analyst should not 
jump too hastily to conclusions that may be valid only for the 
limited number of situations studied ; they could be nullified by 
factors not considered in the particular study. A system study is 
undertaken principally because conclusions are required on which 
to base design of the system. Consequently, the analyst must make 
a particular effort to reach a set of conclusions which are techni- 
cally sound. Then he must present them in such a form that their 
significance will be understood easily by those directing the project. 
Presentation of raw data is usually of very little value as compared 
with presentation of data for a few select cases. They should be 
supplemented by curves and tables to summarize the performance 
of the system under a variety of operating conditions to reflect 
changes in the important system parameters. Unless the analyst 
succeeds in reducing his results to a form readily understandable 
by others, he has not succeeded. The results of system studies are 
usually presented in the form of a report rather than as a working 
piece of equipment. Therefore, the analyst must also master the 
ability to express himself clearly, concisely, and convincingly on 
paper. Although the engineer frequently shies away from the task 
of writing reports, he must learn to crystallize his ideas if he is to 
prepare a report that effectively serves its purpose. 
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Much work remains to be done after the system analyst presents 
his recommendations on the design of a system and before a working 
prototype of the operational system becomes available. If the 
system is at all complicated, the stages outlined in this discussion 
may have to be repeated before all aspects of the system have been 
examined in sufficient detail. In many instances, testing of portions 
of the system is possible long before the complete system design 
is finished. For example, during the development of an anti-air- 
craft guided missile, preliminary systems tests might be made on a 
vehicle which includes only the airframe and a power plant, but 
no provision for controlling the missile. Later, an autopilot might 
be included. Now tests can be made to determine the ability of 
the autopilot to accomplish simple maneuvers which might be 
programmed into the control system. Still later, a vehicle including 
guidance equipment might be tested against a drone target plane. 
But even at this stage in the development, modifications in airframe 
or power plant may be made before a truly operational missile is 
produced. 


MACHINE METHODS 


A wealth of mathematical background exists from which the 
analyst can draw in his studies of the mathematical model used 
to describe a system. Nevertheless, the equations encountered in the 
majority of the really challenging systems problems cannot be 
solved with. presently available analytic techniques. In many 
instances, such equations cannot even be formulated. Consequently, 
the analyst is forced to employ machine methods of computation 
in his study of the equations. Occasionally, he may have to skip 
this phase of the analysis and substitute in its place the study of a 
physical model of the system or experimentation on a full-scale 
system. These latter alternatives are expensive. Also, they do not 
permit as wide a range of parameters to be investigated as is desi- 
rable in a systems study. Fortunately, computational machine 
methods have been developed concurrently with control-systems 
~ techniques. Therefore, the analyst who uses both analytic and 
machine techniques in their proper balance is generally able to 
give a relatively complete description of the performance of the 
proposed systems. The goal of simulation is the effective utilization 
of machine computation methods in order to establish adequate 
criteria for the design and development of proposed systems. 

Analog techniques are well adapted to studies where the analog 
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has a counterpart in the physical system. Aircraft autopilots or 
automobile suspension systems are typical of this situation. When 
used to study such systems, the accuracy of a properly operating 
computer is probably better than the accuracy with which the 
physical system can be constructed and or with which its parame- 
ters can be measured experimentally. In problems where portions 
of the computer have no physical counterpart to reality, the analog 
computer may be hard pressed to yield results with the required 
accuracy. For example, if terrain features significantly in influen- 
cing propagation, the geographical relationships between battle 
field locations are relatively difficult to assess accurately on an 
analog computer but correct solution of these geometrical rela- 
tions is inherent in the real physical situation. 

Most analog computer studies of complex physical systems are 
made on computers operating with approximately a one-to-one 
ratio between the time required for the event to take place in the 
physical system and the time required for the computer to obtain a 
solution. If extremely fast or slow events are to be studied, this 
ratio is sometimes adjusted so that the solution will require appro- 
ximately one minute of machine time, whether the actual event 
happens in ten seconds or ten minutes. With solution times from 
one to two minutes, drift in modern analog equipment introduces 
negligible inaccuracies, and electromechanical computing elements 
and output units such as strip recorders and plotting tables can be 
employed. Sometimes components of the physical system or even 
human operators are included in the computer setup. In such 
situations meaningful results are obtained only if the machine can 
be operated satisfactorily on a one-to-one time scale. In many 
studies of complex systems, relatively little is gained if the time for 
a solution is reduced, for example, from one minute to fifteen 
seconds. The majority of the time spent on an analog computer 
program is used for setup, to check the computer, and to decide 
on necessary changes. After all, the type of results expected from 
the computer or from the modified system must reflect their per- 
formance in the chosen environment. 

Digital computers are employed in the study of mathematical 
systems models. At the present time, the principal technical defi- 
ciency associated with their use in the study of complex systems is 
their speed. An analog machine operates simultaneously on every 
term in a mathematical model of a system. Thus it requires no 
longer to solve a very complex problem than to solve a simple one, 
although it is virtually impossible to assign an accuracy figure to 
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some of the computations. On the other hand, the time required 
to obtain a solution on a digital machine increases with the com- 
plexity of the problem. Even with the large, high-speed, electronic 
digital machines now available, the ratio of solution time to real 
time may run high for moderately complex systems problems. 

Nevertheless, digital machines at present play an important 
rôle when only a few solutions are required or when accuracy requi- 
rements are high. For example, digital solutions for one or two 
cases from a complete family of solutions may prove extremely 
valuable as check points in an extensive analog study. Further- 
more, after completion of an analog study, the analyst may find 
that his over-all understanding of the system would be greater if a 
few additional case studies were available. In this instance, it may 
be more efficient to obtain these few solutions digitally, particu- 
larly if the basic digital program had been coded in the process 
of obtaining a digital check solution prior to the analog work. 

The speed of digital machines has been increased significantiy 
in the past decade. It is quite possible that in another ten years the 
speed advantage now held by analog machines in the solution of 
complex sets of differential equations may be reduced to a point 
where, on this basis alone, analog machines will have lost their 
present advantage. However, as the two types of machines are now 
operated, the engineer generally feels more at home with an analog 
machine ; a digital machine appeals more to a mathematician. The 
engineer is accustomed to working with block diagrams. The very 
close analogy between the block diagram that he uses to describe 
the system under investigation and the computer setup diagram 
offers distinct advantages. If the engineer who is analyzing the 
system follows the study through its various stages and takes actual 
responsibility for operating the analog computer, then he can apply 
all his knowledge of the system’s operation to make on-the-spot 
decisions concerning the desired parameter values. From his first- 
hand observations of difficulties with the analog setup, he may be 
able to predict weaknesses that would arise in the analogous phy- 
sical system. 

This close association between the systems engineer and his com- 
puting machine does not always exist when a digital machine is 
used. Usually, the system analyst prepares a set of equations which 
he gives to a programmer who prepares a program for the operation 
of the digital computer. The programmer in turn gives the program 
to someone who punches a set of cards, or a paper tape, by which 
the problem is fed into the computer. Still a fourth person handles 
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actual operation of the machine. In most cases, the machine opera- 
tor has little knowledge of, or interest in, the specific problem on 
the machine. On the basis of the efficiency with which each step 
can be performed, this arrangement is ideal, but it is apt to be 
inefficient in the over-all sense of providing maximum assistance 
to the analyst in the design of his system. 


SIMULATION METHODOLOGY 


The formulation of mathematical models is an intuitive process 
and their validation is not always feasible or even desirable, as in 
war gaming. Turning now to the problem of model implementation 
(for example, on digital computers) we shall look at some of the 
techniques that have been developed for this purpose. 

The first section of this chapter will deal with analytical techni- 
ques capable of yielding models of any desired precision. As is the 
case with numerical differentiation or integration (Runge-Kutta 
method), the finer the granulation, the more accurate will be the 
results obtained from the process. However, the amount of work to 
be done will also increase. Relevant analytical techniques discussed 
will include three distinct levels of time granulation ; micro, sem1- 
macro, and macro-simulation. 

Considerable work has also been done in the area of programming. 
The earliest programs required that simulated events be distinct 
and sequential (Sequential Processing = SP). However, two very 
powerful techniques came soon into being and have been used 
extensively where event simultaneity is a factor. They are some- 
times referred to as Status Vector Cycling (SVC) and Time Status 
Recording (TSR). Synonymous terms are synchronous and asyn- 
chronous (‘* most imminent event technique ’’). 

The use of simulation programs has also been studied extensively. 
Controlled experiments with randomly chosen events in a stochastic 
process are known as Monte Carlo experiments. Emphasis is now 
being placed on the concept of statistical control of the random 
numbers which set the pace of the program. Therefore, a brief 
mention will be made of the so-called Monte Carlo swindle; also, of 
the problem of when to terminate a single run or a Monte Carlo 
experiment. 


ANALYTICAL TECHNIQUES 


We may distinguish three levels in the digital simulation of 
(electronic) systems. These levels are characterized by the fineness 
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of time granulation and the degree to which individual system 
functions are simulated and analyzed. 


Micro-simulation. 


This term implies that the time quanta used are extremely small ; 
how small depends upon the desired degree of precision of which the 
simulation is to be capable. The action of single electronic compo- 
nents may be simulated by separate computer subroutines. Elec- 
tronic stress of individual components is then assessed statistically. 
As a result, we may obtain, for example, the mean time between 
failures (MTBF) for the system by aggregation of MTBF's for single 
components. Present component densities of several thousand 
units per cubic foot limit the scope of micro-simulation to at most 
a few drawers, or racks, of equipment. Time quanta must be chosen 
small enough to permit evaluation of pulse rise time, noise, and 
pulse decay in the simulated digital equipment. Micro-simulation 
permits the evaluation of equipment reliability, failure rates, signal 
response, or of the probability of inadvertent generation of pulses 
or signals. 


Semi-macro-simulation. 


Going to a higher order of time aggregation, this form of simu- 
lation analyzes pulses and cycles in a system. The applicable time 
quanta in synchronous or asynchronous systems are chosen to equal 
the basic time of one pulse or cycle. The simulation program may 
provide for a mapping of all memory bits in the system such as 
shift registers, flipflops, gates, cores, etc. Error rates may be assi- 
gned to individual bit functions. They are not derived from a mi- 
croscopic simulation but are supplied in parametric form. For 
example, the error rate applying to a shift register may be given 
as the probability for inverting a mark or a space within each ele- 
ment. 

Semi-macro-simulation is limited in scope to medium size systems, 
such as digital computers, message switching devices, or the like. 
If individual error rates are available, this type of simulation may 
be used in two fundamentally different ways. First, we may check 
consistency and completeness of the logical rules which were incor- 
porated into the design of the equipment : operation without errors. 
Secondly, we may determine. system operational characteristics 
in the presence of given or postulated error rates (or noise). 
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Macro-simulation. 


This type of simulation is again less finely granulated than the 
semi-macro-simulation. Dependent upon the form chosen, the 
granulation may be by messages, or by time blocks. The method 
applies to large-scale communication networks, for example. It 
provides a tool to validate a given network design by comparing 
message queues, call rates, lost calls, etc., with requirements expres- 
sed in form of dynamic communication needs or system constraints. 
It also provides a means of evaluating the system’s response to 
local or over-all transients. ‘* Local’’ transients might reflect 
inputs for periods of pre-hostility or war, in cases of area-limited 
conflicts ; over-all transients would occur in the global conflict case. 

System operations in the presence of errors (or noise) may also 
be incorporated into such a program. If the time quantum is chosen 
on a message-by-message basis, degradation of digital or voice 
traffic from link to link can be simulated. Of course, one must 
postulate separate models for the compounding of this degradation 
if traffic is carried over more than one link. However, such models 
are now available both for digital and for voice traffic. 

Macro-simulation may be carried out either on a one-to-one basis 
or on a Statistical basis. In the one-to-one simulation the identity 
of individual calls, or digital messages, is retained and recorded. 
The combination of this activity into call rates leads to statistical 
simulation, where the identity of individual calls is neither observed 
nor recorded. Experts have contended for some time that it is extre- 
mely difficult, especially in the case of saturated systems, to vali- 
date the transition from the one-to-one to the statistical model. 
Indeed, great conceptual difficulties arise in connection with this 
transition. The statistical parameters put to work can be at best 
empirical estimates obtained either from field experience or from a 
one-to-one simulation model. 


PROGRAMMING TECHNIQUES 


Three different techniques which have been used successfully in 
the simulation of large systems will now be discussed. 


Sequential processing (SP). 


There exist systems in which no more than one activity can be 
initiated at any one time while the duration of each activity is not 
necessarily always the same. A typical system of such a nature is a 
sequential, digital computer. 
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Simulation by sequential processing assigns to every activity 
a time of initiation and a time of termination. No activity may 
be initiated unless the preceding activity has been terminated. 
The choice of an activity depends upon the system state existing 
at the termination of the preceding activity. For example, in simu- 
lating a digital computer, execution of one instruction may cause 
a decision as to which instruction is to be carried out next. This 
would be the case for transfer and compare instructions simulated 
by the program. 

Simulation by sequential processing is limited in scope to those 
devices and systems which act similarly to a sequential, digital 
computer. A communications network permits simultaneous injec- 
tion of more than one message, at more than one point of the system. 
Therefore, it does not admit of sequential processing in the strict 
sense. Yet, certain subcomponents in a communication network 
could be simulated realistically by sequential processing while the 
overall system permits simultaneity of activities in several areas. 


Status vector cycling (SVC). 


Most large systems may be defined conceptually by a time- 
dependent status vector. The elements of this vector include the 
state of all digital and analog devices in the system, at any given 
time. Analog devices may be capable of assuming more than two 
states. They would take their assigned values from a continuous 
function (which in turn represents the model for the device). Bista- 
ble components may be assigned values zero or one (mutually 
exclusive) at any given time. 

Holding the state of all devices and components of the system 
in suspense, for a short moment of time, yields the system state 
vector, for that instant of time. A short time later, some or all 
components of this vector may have changed their status. Thus a 
short time later another state vector obtains. The system state 
vectors form a sequence in time which reflects the reaction of the 
system to external inputs and internal feedback. 

If the cycle time for setting up this instantaneous system status 

vector is chosen small enough, only a few of the vector elements 
will change from one cycle to the next. The nature of these changes 
may be deterministic or stochastic. A typical example for deter- 
ministic changes is afforded by the actions of a mechanical relay, of 
a binary half-adder, or the like. Given the state of the device during 
one cycle and certain inputs, the next state of the device is uni- 
quely determined. An example of probabilistic element changes in 
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the state vector is afforded by the termination of a voice call in a 
telephone network. In this case the probability that a line which 
is busy will become free during the next system cycle has some 
preassigned (Bernoullian) probability which is applied to the call 
at each cycle, as long as that particular line is active. In effect, it 
creates binomially distributed random holding times which ap- 
proach closely exponential holding times if the binomial proba- 
bility is chosen small enough. 

For small time granulations few changes occur in the elements of 
the system state vector. However, during each cycle all elements 
capable of probabilistic state changes must be operated upon. The 
control program logic for such programs is very simple, but an 
enormous amount of work is performed during the computer simu- 
lation. For example, if the system cycle were three seconds and 
the average call length were two minutes, on the average forty 
probabilistic calculations are made between initiation and termi- 
nation of one such call. 


Time status recording (TSR). 


The most efficient and economical method of simulating systems 
in which more than one event may take place simultaneously is 
afforded by the time-honored TSR technique. Instead of perfor- 
ming a large number of Bernoulli-trials to determine initiation or 
termination of a single event, the times of initiation and termination 
only are determined from applicable distributions (Poisson, expo- 
nential, etc.). Deterministic reactions caused by the event itself 
are carried forward, in time, to all other logically connected events. 
Thus, with every event there is associated an action time. This 
action time, as well as the next action time for the same event, is 
recorded in a time status vector. 

For example, if voice calls are initiated at some given rate, their 
time of initiation is recorded. Simultaneously their time of termi- 
nation is noted and is followed by another call initiation time, etc. 
The fundamental distinction between TSR and SVC programs, 
therefore, lies in the number of times probabilistic subroutines are 
addressed. 

The distinctive feature of the TSR control program is the sorting 
of events by action time. A fictitious clock is set up in the computer 
to determine which event must be processed next. Upon processing 
of this next event its deterministic influence on other events is 
noted in all logically connected time status subrecords. The com- 
puter program holds a tight rein on the listing of all system ele- 
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ments which are subject to seizure, release, or analog reactions. 
It does this within the framework of a large number of time status 
subrecords. 

It is readily seen that the TSR method of simulation occupies no 
more memory space for the mapping of system functions and ele- 
ments are addressed once, and only once, upon a change of status. 
In addition, a time sort follows the completion of each action to 
ensure that the next element processed is indeed that which comes 
next in time. 


MONTE CARLO EXPERIMENTS 


Digital computer simulation programs become increasingly 
dependent upon Monte Carlo experiments to determine simulta- 
neously a multiplicity of systems parameters. By empirical methods, 
the design engineer may thus obtain parameter values which could 
not be found by analytic solution. In many instances, mathema- 
tical formulation of the basic problems is not even feasible. 

In the course of his work, the design engineer may set up an 
acceptable model for his system and subject it to computer simu- 
lation. The results obtained from this simulation do not really apply 
to the system ; they apply only to the model chosen. The model 
must be validated against the system, by field tests or selected expe- 
rimentation. Once validated, the model can be accepted as a repre- 
sentation of the real system ; it is then used, by extension, to obtain 
information which exceeds the scope of field testing. 

Two fundamental problems arise prior or during the running 
of such simulation models on digital computers : (1) when to stop 
a single run, (ii) when to stop a Monte Carlo (multi-run) experi- 
ment. The first problem concerns the detection of transients in the 
simulated system and the time at which the system may assumed 
to be in steady state. The second problem concerns the reliability 
that may be attached to the estimated system parameter(s). 

We observe, in passing, that system simulation does not cause 
a run termination problem if the simulated, real time span has 
well-defined limits. This case may occur in the simulation of missile 
and aircraft flights, where the simulated real time extends from a 
launch-phase, through mid-course guidance (if any), to a terminal 
phase involving impact or landing predictions. 

The general approach to the Monte Carlo analysis should be: 
how much information is available for each item in the state 
vector ? The decision to terminate a single run, or a Monte Carlo 
experiment, must then take into account the quantity and quality 
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of this information. Frequently, it is not available before running 
the computer program. As a matter of fact, if all of this information 
were available, there would be little need to run the simulation. 
Detailed knowledge about the initial transients, especially, can 
hardly become available without running the program. Therefore, 
one reason for the running of such programs is the need to deter- 
mine the duration of the transient and to fix a time at which steady 
state of the system may be assumed. 


The run termination problem. 


The program is only a system model of assumed validity. Its 
operation on the digital computer, on the other hand, is a well- 
defined pseudo-stochastic process. While the engineer looks at it 
as a simulation of his system, the program analyst views it as an 
input-output problem. He is intimately familiar with the design of 
the program. He “ knows ” that each run will converge towards 
some steady state, despite initial transients. Of course, the ‘‘ steady 
state ’’ may itself be a stochastic process with variations or fluctua- 
tions on certain elements in the state vector. But as viewed by the 
program analyst it is a predictable process, even though the ele- 
ments of the state vector may exhibit limited ‘‘ excursions ”’. 

In view of the above, the following two cases must be considered : 
(i) assumption of no initial transient, (ii) definite knowledge, or, 
assumption of existence of initial transient. In the first case, a 
statistical test must be applied to the ‘‘ weakest ’’ element of the 
state vector. This test determines whether or not there was an 
initial transient. In the second case, the same test may be used 
to estimate the termination time of the (initial) transient and to 
establish some time beyond which the system may be viewed as 
being in steady state. | 

Monte Carlo experiments are used to study the simulated system 
by making repeated runs, on the basis of an experimental design. 
Input parameters are varied under the control of randomized expe- 
riments ; random numbers are used to control stochastic processes. 
These two elements, then, determine the over-all variance for the 
program outputs. 

If only a single parameter had to be analyzed by simulation, 
certain advance estimates could be made with regard to the number 
of runs needed to obtain an estimate for this parameter with prede- 
termined precision and confidence. However, the situation is usually 
not quite as simple as all that. First, the use of a large-scale simu- 
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lation program for the purpose of estimating a single systems para- 
meter would be rather inefficient. The primary reason for writing 
such programs, in the first place, is the need for simultaneous, 
multiple parameter estimation. A single parameter would likely 
require a much simpler program. A more sophisticated reason for 
the problem complexity lies in the fact that there usually does not 


exist any prior knowledge about the statistical distribution of the 
parameters. 


Statistical control. 


We may wish to exercise statistical control over a Monte Carlo 
experiment by terminating it as soon as a certain parameter esti- 
mate can be made with prescribed confidence and precision. A 
repetitive procedure could be developed for that purpose. This 
procedure monitors the course of the Monte Carlo experiment, like 
an invisible observer. It makes continual re-estimates of the para- 
meter, after each run within the Monte Carlo experiment. When 
these re-estimates indicate that the desired precision and confidence 
have been attained, the control procedure causes the program to 
stop. 

Such a Monte Carlo control program is illustrated by the flow 
chart in Fig. 1. Following program initialization, the number of 
runs made thus far is counted and entered into storage. The accu- 
mulated data are then tested for normality. If they don’t pass this 
test, the observer is left to his own devices as to how he should 
proceed. Recommended procedure might be to make another five 
to ten runs and to try again. If the data do pass the normality test, 
an estimate N is obtained for the total number of runs required to 
attain a parameter estimate of given precision and confidence. 

At this point, the analyst may have selected one of two options. 
Either he wishes to observe the Monte Carlo experiment on a run- 
by-run basis, or he wants to get on with the business and finish the 
program as fast as possible. In the former case, he may even use 
the cathode-ray tube to display the run-by-run progress. In the 
latter’ case, another N-n runs are made, and the parameter esti- 
mate is only then tested and repeated runs are made, if needed. 


The Monte Carlo swindle. 


One type of manipulation that is commonly used to decrease any 
error variance is known as the Monte Carlo swindle. As an 
example, consider the use of the original random number inputs. 
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After a design has been constructed, each run could be started 
with the same set of random numbers: the perfect Monte Carlo 
swindle! By this method, an artificial correlation is generated 
between runs, suppressing much of the error variance. 

Another type of manipulation that has been used to decrease the 
error variance is the technique of sampling without replacement 
from the universe of random numbers. In this case, the expected 
distribution of events is predetermined and is placed in the model 
before the run; drawings without replacement are then made 
randomly. Hence, any variance which accrues depends largely 
on the order in which the events happen. 


NEXT SIMULATION RUN 
WITHIN THE MONTE 
CARLO EXPERIMENT 


START WITH INITIAL 
NUMBER OF MONTE 
CARLO EXPERIMENTS 


MAKE OTHER 
DECISION ON 
NUMBER OF RUN 
NEEDED 


yes 


ESTIMATE N, THE 
TOTAL NUMBER OF 
RUNS REQUIRED 
TO ATTAIN DESIRED 
PRECISION AND 
CONFIDENCE 


Fic. 1.— Monte Carlo experiment control program. 


It is evident that the special characteristics of every simulation 
model require the careful study of experimental designs. Until more 
theoretical work has been done, investigators must be content with 
experimentation in model design, and with the acute problem of 
accomplishing this by the shortest method, in the least amount of 
time, at minimum cost. 
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EXAMPLES OF SIMULATION PROGRAMS 


The thesis has been proposed that all the work done on digital 
computers is merely simulation — and indeed it is! Man has the 
capacity to grasp and manipulate the concept of infinity, but he 
cannot possibly attain it through finite processes. This thought now 
suggests the more mundane discussion of where and when digital 
simulation programs have been used. 

One area in which digital simulation has achieved notable success 
concerns the design and development of engineering devices. Cer- 
tain simulation programs (for example, UDOFT, the Universal 
Digital Operational Flight Trainer) require special purpose com- 
puters ; in turn these computers may be simulated on other com- 
puters, prior to their construction (example : PSEUDOFT on IBM 
704) (1). Other programs are written for general purpose computers 
such as most programs written for nuclear reactor simulation. In 
fact, this class of programs causes an interesting observation. It is 
imperative that one has a good understanding of the behavior of 
such systems in advance : if the ‘* plumbing ’’, so to speak, would 
require ‘* fixing ’’ after the construction of the reactor, we would 
soon run out of plumbers willing to sacrifice their lives. Other engi- 
neering devices which are extensively simulated at great savings 
to the taxpayer are missiles. Since the price per shot may run as 
high as many million dollars, one can imagine the cost of these 
development programs without simulation, since it is not unusual 
to ‘‘ fire ’’ thousands of ICBM’s in the course of one computer run. 

A second area where simulation has been instrumental in provi- 
ding excellent services is the simulation of one computer on another. 
These programs may use sequential processing since the basic part 
of the computer logic is usually nonsimultaneous in nature. They 
are written for two reasons. First, design engineers may wish to 
check out the logic of a new computer, prior to its construction ; or, 
programmers may wish to debug subroutines before the hardware 

becomes available. This type of computer simulation might be 
considered as internal, since it involves only noncompetitive aspects. 
A second, more compelling reason for writing computer simulation 
programs has competitive aspects. As new computers with very 
fast processing speeds become available, it may be desirable to 


(1) Relevant work was done at the Waltham Laboratories of Sylvania Electronic 
Systems. 
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avoid reprogramming and to use old programs on the new machine. 
In such cases, one must meet two objectives : (i) the simulation 
program must be at least as fast as the old (object) computer, and 
(ii) it must run cheaper on the new (subject) computer (say, on a 
cost-per-bit or cost-per-unit-of-time basis). A classical example 
is the simulation of IBM 650 on IBM 704. 

Simulation programs have been extremely useful to the commu- 
nications industry. The designers of communications networks, 
especially of global systems, could not possibly build many of them 
and then discard all but a few, e. g., keep the ‘* good ”’ systems. 
Such a procedure would soon exhaust the world’s resources. The- 
refore, communications system engineers have come to rely on 
simulation as their most important tool. Appropriate techniques 
and pertinent model types will differ from one case to the next, 
especially where voice or digital nets are concerned. For example, 
consider global nets for the military (Aircom, Globecom, Unicom, 
Navcom, etc.) and commercial networks (AT&T, RCA, etc.), all 
planned or operational ! 

Finally, there is one area which has furnished more food for 
thought and simulation programs than any other: game theory. 
The term originated with Oskar Morgenstern and the late John von 
Neumann who wrote the classic on the subject [9]. More recently, 
the term has become almost synonymous with war games. It assu- 
mes, in people’s minds at least, all the sinister aspects of large 
weapons systems and mass destruction. Typical war games are 
being simulated on various scales. Small battles are being ‘* fought ”” 
with conventional weapons ; global battles are ‘‘ waged ’”’ with 
planes carrying nuclear bombs, with ICBM’s being launched both 
singly and in salvos. The damage to the various parties is assessed 
by statistical summary programs at the conclusion of individual 
runs. One of the better known large games of this type is the Air 
Battle Model (ABM) which was programmed by Technical Opera- 
tions Inc. and the subsequent Damage Assessment Model (ADAM). 

We must not leave this topic without mentioning that gaming 
is not always war gaming. In fact, computer programs are being 
used increasingly for nondestructive purposes, by playing ‘* peace 
games ’’. In our competitive society game theory provides a very 
useful tool to management in the area of decision-making on a 
grand scale. Management games are frequently played for educa- 
tional purposes only. However, several operational programs have 
been written to furnish control of large contracts and procurements, 
to pinpoint possible trouble spots in advance, and to find where 
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delays could cause more or less serious difficulties in the over-all 
program. PERT (Program Evaluation Review Technique) is such 
a management technique which assists with the reduction of re- 
search and development time on large systems contracts. 

Thus we see that simulation programs have become a vital tool 
in diversified areas ranging from straight engineering to man-ma- 
chine systems. The rôle played by the computer is that of a control 
device. Peace games, in particular, are making increasingly impor- 
tant contributions to man’s better understanding of his environ- 
ment. They furnish him also a better control of factors which he 
has never before been able to understand. 


COMMUNICATIONS NETWORKS 


One of the major purposes of any proposed simulation program 
is the evaluation of existing or the validation of proposed system 
designs. The program may also have to assess the incremental 
utility of equipment changes or of system augmentations. There- 
fore, emphasis must be placed upon two parts supplementary to 
the program : first, the initialization phase, which fixes input 
parameters, operating conditions, and applicable war gaming ; 
second, the post-processing phase, which reduces to a readily under- 
standable format the mass of data accruing from the simulation 
program. This part of the program, in particular, makes extensive 
use of the concept of ‘ management by exception ”’. It should not 
produce reams of printed outputs of a consistently stagnant or 
noisy character, but it should exhibit concisely only those output 
items of essential information character. 

Single runs with a communication network simulator are usually 
made to obtain information on selected system parameters which 
have the nature of random variables. Grossly, they may be classified 
as follows : 


(i). Integral valued random variables, usually non-negative. 
Examples are message counts, missile counts, queue sizes, lost call 
counts, etc. 

(ii) Continuous random variables of a primary nature, usually 
(but not always) non-negative. They may include single call delays, 
miss distances of missiles, arrival times of messages, etc. 

(iii) Continuous, dimensional random variables of a derived 
nature. They may include average call delays, average queue sizes 
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(a queue is then assumed to have ‘* length ’’), mean time between 
failures, etc. 

(iv) Continuous, nondimensional random variables of a derived 
nature. Examples are relative frequencies, probabilities, fraction 
of a missile force delivered on target, failure rates, lost call rates, etc. 


The. broad purpose of communication system simulation includes 
usually the evaluation of more than one system parameter. There- 
fore, most single runs are made with the aim to produce a sequence 
of state vectors at wnich the analyst ‘* looks ’”’ to determine tran- 
sients or steady state. The entries in these state vectors may be any 
of the aforementioned random variables. During a single run, some 
of the entries may reach steady state earlier than others. This is 
particularly true when simulating a communication network, not 
wholly saturated with traffic. If traffic rates in parts of the system 
differ appreciably from one another, some system parts with lower 
rates may attain steady state earlier than others. 

In the following subsections we shall describe in some detail the 
three phases of such programs : initialization, simulation, and post- 
processing. 


INITIALIZATION PHASE 


This part of the program must establish and /or set all the param- 
eters required for the subsequent simulation run. It supplies the 
simulation program with network connectivity matrix, applicable 
error rates, times of singular events, etc. For example, in the case 
of dynamic war gaming, the simulated attack may consist of several 
waves arriving at different times and/or at different locations. 

The modular format of the simulation program must permit its 
organic growth exactly as it would occur in the real, global commu- 
nications system. Thus, the connectivity matrix must be capable of 
expressing line and trunk patches, or of permanently allocated 
channels. It must facilitate simulation of manual, semi-automatic, 
and automatic switching. The simplest way in which this organic 
growth of connectivity and switching can be accomplished is to 
consider the end point of any communications channel or line as an 
interface. This interface belongs to a certain item in the time- 
status-record and is permanently connected to some other inter- 
face ; or it is connected to one or more other interfaces by switching 
or patching operations. 

War gaming considerations require separate identification of 
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individual communication links, channel groups, or trunks. The 
earlier introduced concept of interfaces at the end points of each 
line, trunk, or channel serves to identify individual communication 
links and permits the assignment of a reliability figure to each ele- 
ment. Thus, in the case of jamming or other covert enemy action, 
this reliability figure may be adjusted according to some computer- 
stored model. 

Another item which must be given thought during the initializ- 
ation phase is the implementation of different routing doctrines. 


INITIALIZATION PHASE : 
1.SET INPUT OPTIONS 
2.SET PROGRAM SWITCHES 
3.SET ON-LINE STATISTICS 
4.SET OUTPUT OPTIONS 

5. START TIMER 


FIND “NEXT” ITEM IN TSR 
TIME CONTROL LISTING 


SELECT “FIRST ITEM IN TSR 
CONTROL LISTING 


START PROCESSING OF THIS 

ACTION ITEM IN TSR: 

1. COLLECT INPUTS RECEIVED 
FROM EARLIER ACTION ITEMS, 
IF ANY. 

2.PROCESS THIS ACTION ITEM 

3.APPLY ERRORS, DETERIORATION 
RATES, IF ANY. 

4.DISTRIBUTE OUTPUTS AND 
“ARRIVAL” TIMES. 

5.OBTAIN SUCCEEDING (NEW) 

ACTION TIME FOR THIS ITEM. 


ENTER THIS NEW ACTION ITEM 
WITH ITS NEW ACTION TIME 
INTO TSR CONTROL LISTING; 
DELETE OLD ITEM 


UPDATE ON-LINE STATISTICS; 
UPDATE DISPLAY, IF ANY 


WRITE OUTPUT TAPE ITEMS 


POST-PROCESSING PHASE: 

1. PREPARE OUTPUT TAPE FOR 
PRINTOUT 

2.CONDENSE OUTPUT DATA 

STATISTICALLY 


END OF REAL-TIME RUN? 


Fic. 2. — Overall flowchart for TSR macro-simulation. 
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The simulation program may have to permit assessment of such 
radically different routing techniques as “‘ direct distance dialing 3 
(DDD), or ‘* store-and-forward ”’. 

The introduction of error rates into the simulation program must 
also take place during initialization. These error rates are of different 
types. Voice traffic over several links in tandem suffers essentially 
additive degradation. It is based upon the characteristics of the 
links over which the traffic is routed. Digitalized traffic, on the 
other hand, even if it is transmitted over a nominal 4 kc channel 
is subject to error rates which cause the inversion of single message 
bits. In macro-simulation, sometimes only cumulative effects can be 
taken into account. A one-to-one scheme must also retain the iden- 
tity of individual messages. In that case, it must be possible to 
simulate the receipt of a message at any intermediate or at the 
ultimate location and to assign to the message output a compound 
degradation coefficient. Voice messages would then be considered 
unintelligible if this degradation coefficient exceeds a certain thresh- 
old. For digital traffic, various thresholds for message intelligi- 
bility will be required to cope with teletype, speech compression, 
data coding, crypto, etc. 


SIMULATION PHASE 


Fig. 2 illustrates the over-all flow chart for the processing of 
time-status-records. The principal feature of this technique is the 
processing of ‘* action items ”’ at the time of their activation. This 
‘action time ’’ in an action item (‘* most imminent event ’’) 
implies the immediate and the next action time. The action items 
are stored in a large table ; they are sorted by action time. Once 
an item has been acted upon it is removed from the table and the 
next item in order of time is processed. Removal of an item from 
this list, however, may cause the item to reappear somewhere else 
in the list, for later action. For example, for a voice call as a typical 
action item, its first action time coincides with the time of its ini- 
tiation ; its second action time with the time of its termination. 

Processing of an action item at its action time may cause deletion 
of other action items, changes in their action times, or the creation 
of new action items with their associated action times. For example, 
consider again the case of a voice call : it is actually composed of a 
large number of sequential action items. In a semi-automatic 
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system, as shown in Fig. 3, the following events take place prior to 
the completion of a connection between two subscribers : 


(a) Subscriber’s call initiation. 
(b) Initial local switchboard response. 
(c) Initial local line selection. 
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Fic. 3. — Typical, semi-automatic message switching network. 
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Trunk switchboard response. 

First trunk channel selection. 
Automatic message switching. 
Second trunk channel selection. 
Automatic message switching. 

Third trunk channel selection. 
Automatic message switching. 
Terminal trunk switchboard response. 
Terminal local switchboard response. 
Terminal local line selection. 
Terminal subscriber unit response. 
Terminal subscriber response. 
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In passing, we note that the holding time for a voice call cannot be 
applied to the first action item listed above until the circuit has 
been completed and the second subscriber has answered the call. 
Thus, this turns out to be a typical example of the modification of 
the action time of one item by the processing of another 

The table of action items is composed of two types of entries. 
First, there are permanent entries for such items as call or message 
generating sources and sinks. This part of the table is exactly as 
large as the number of subscribers. Prior to initialization, all of 
these subscribers are ‘ not busy ’’. Another part of the table con- 
tains generated items. The latter belong to the class of chain events 
which are created by ‘‘ busy ”’ action items of the first class. Typical 
chain events include activation and de-activation of dial tones, 
lines, switching devices, etc. 

It is trivial that in a network completely saturated with traffic 
the number of action items in the list depends on the number of 
action items generated and on the maximum number of busy 
customers of which the saturated network permits. If the system is 
overdesigned, all customers may be busy at once ; if it is not over- 
designed, queueing is forced on some of the customers. The total 
listing of action items is thus considerably larger than the number 
of customers, in either case. 


POST-PROCESSING PHASE 


Certain statistical information will be accumulated “ on-line ’’, 
while running the simulation program. The flow chart, Fig. 2, also 
indicates the possibility of providing on-line information to output 
devices such as a cathode-ray display tube. For the sake of keeping 
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the program running time low, on-line processing should be held 
to a minimum. 

The bulk of the data processing work load should always be 
carried in the post-processing program. Here, data deposited on 
magnetic tape files during the simulation run may be scanned and 
processed without interfering with the real time aspects of the 
simulation program. The post-processing phase must permit the 
evaluation of any desired effects, such as choice of routing doctri- 
nes, static or dynamic war gaming, network survivability and 
vulnerability, etc. In particular, the post-processor must reduce 
this statistical information into indices to validate system designs 
and to evaluate the incremental utility derived from any proposed 
network augmentations. 

The construction of such indices is a nontrivial statistical pro- 
blem. Global communication systems are not characterized by 
‘* average number of lost calls ’’, ‘ average message delay ”’, etc. 
These averages lose their meaning as the networks become larger. 
For example, if one-half of the network were completely saturated 
and the other half completely quiescent, the ‘‘ average ’’ load 
figure for the system would be 50 per cent. It is easily seen that such 
a figure has no real meaning. Instead, traffic load matrices should 
be set up, on a link-by-link or on a regional basis. Comparative 
evaluation is then furnished by statistical comparisons, rather than 
single indices. 


COMMUNICATIONS NETWORKS SIMULATION PROGRAM 


We shall now demonstrate how to simulate a communications 
network on a digital computer, using the method of status vector 
recording. Naturally, it will be a relatively simple example. In 
fact, we shall restrict ourselves to an example which does not require 
the use of a digital computer ; nevertheless, it will reflect most of 
the properties of a real-life simulation program. In other words, 
what we are about to do is to simulate a simulation program. 

We shall first describe the physical system, its operational cha- 
racteristics, and the chosen routing doctrine. Next, we develop a 
simulation flow chart, paying particular attention to the design 
and contents of the state vector. To implement random operation 
where required, we shall examine the use of random numbers, 
especially if we are to avoid such traps as the earlier cited Monte 
Carlo swindle. In the actual process of simulation, we will then 
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assess the system state by counting lost calls and other pertinent 
variables. Finally, we shall summarize the results in graphical form, 
to capture transients and steady state of the system. 


SYSTEM DESCRIPTION 


Fig. 4 illustrates a simple switching network. Four automatic 
switching centers are arranged in a square. They are connected to 
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Fic. 4. — Physical system description. 
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one another by trunk lines, each of which is capable of carrying 
three simultaneous conversations. In other words, the automatic 
message switching part of the system has three channels. Each 
switching center services two manually operated switchboards over 
two channels. The switchboards have three subscribers. All told, 
the nodal part of the system comprises four switching centers, 
eight switchboards, and twenty-four subscribers. For the purpose 
of identification, the switching centers are numbered from 1 to 4; 
the switchboards from A to H ; and the subscribers from 1 to 24. 
For example, Fig. 4 indicates that subscribers 7, 8, and 9 are con- 
nected to switchboard C, which in turn is connected to the auto- 
matic switching center 2. 

_ To simplify the later demonstration, we shall assume instanta- 
neous operation of all equipments. Therefore, no delays are incurred 
when a subscriber initiates a call, when his partner answers, or 
when he gets a busy signal. Of course, in a more sophisticated 
program, one would incorporate the delays described in this article 
under the title “ simulation phase ”. 


SYSTEM OPERATION 


We postulate two sets of rules : one for the action and reaction of 
subscribers ; the other for the routing of calls through the system. 
Initiation of calls in a realistic telephone network obeys rather 
complex laws. In this simple demonstration, we shall assume that 
the time interval between successive actions by any subscriber has a 
rectangular, or uniform, time distribution, ranging from 00 to gg units. 
For example, it will appear equally likely that subscriber 21 desires 
to resume his activities after 59 time units as after 77 time units, 
following his preceding action ; these numbers were arbitrarily 
chosen from the set of available random numbers (Table 2). 

Actions by the subscribers may be of two kinds : initiation (or the 
desire to initiate) of new calls ; or, termination of existing calls. If 
_ a subscriber is not busy, he goes into action by. trying to initiate a 
call ; if he is busy, at such an action time, he will then terminate an 
existing call. 

We shall also assume that the choice of partners is random. Thus, 
if subscriber 21 wishes to go into action, we shall toss some sort of 
an hypothetical die with 24 faces to determine randomly which 
partner will be assigned to him. Again, this is an oversimplification. 
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It is known that the propensity to interact, between subscribers, is 
inversely proportional to some power of the distance between 
them [6]. 

Next, we set down the routing doctrine to govern assignment of 
lines, channels, trunks, and switchboards. Fig. 4 indicates that 
there exist at most two alternate routes between subscribers. For 
example, if subscriber 23 wishes to call subscriber 9, then only the 
following two alternate routes are available : 23 H 12 C 9 or, 
23 H 1432 C 9. For some other subscriber pairs, only one route can 
be found : 21 G 4 F 18. In the simulation, the following rules will 


apply : 

(i) The primary route shall be the ‘* shortest ’’ route between 
any two subscribers ; in the case of two alternate routes of ‘* equal ”’ 
length, the primary route will be that one going in clockwise direc- 
tion through the automatic switching centers. 

(ii) If there exists an alternate route, it will be considered as the 
secondary route. 


oD 


An example: For subscribers 23 and 10, the primary route is: 
23 H 123 D'i0; the alternate route is: 23 H 143 D 10. Ti all three 
channels between I and 2 were busy with earlier calls, a call from 
23 to 10 would require alternate routing. But if also all three chan- 
nels between switching centers 3 and 4 were occupied, a call from 23 
to 10 could not be completed ; it would become a lost call. This 
event would then be recorded as such in the system state vector. 


PROGRAM FLOW CHART 


A flow chart, similar to that which would be used to program 
for a digital computer simulation, is shown in Fig. 5. We shall 
briefly discuss its several branches. 

The main branch, boxes I-11, covers the majority of call initia- 
tions and terminations. Its prime element is the sorted time status 
record. Here we record the times at which (future) actions of indi- 
vidual subscribers are expected to take place. The actual entry 
into the record are two-tuples such as (13,7), (18,17), ..., etc., to 
indicate that subscriber 7 must take action at time 13 and that 
no further action take place until time 18 when subscriber 17 acts 
(Table 2). In box 2, we select the next item (‘ most imminent 
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event ’’) in the time status record. In box 3, we test the “‘ age ”’ of 
this future item as we do not want to run the simulation program 
beyond some final value. This box provides the run termination time. 
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Fic. 5. — Flow chart for simulation program. 


If the action item does not cause a program stop, it may be 
either a call initiation or a call termination (box 4). If it is a call 
initiation, we will assign this call a random partner (box 5) and 
select its primary routing (box 6). If this circuit can be completed 
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(box 7), and if the assigned partner is not busy (box 8), the call is 
effectively placed into the system, thereby increasing the system’s 
busyness (box 9). Since this was a call initiation, we must select a 
future (random) time for the termination of this call (box 10). This 
call termination time is entered into the time status record (box IT) 
to complete the main leg of the program. 

The program has several minor legs. Run termination is caused by 
box 3. Thence we go to box 12 and prepare a statistical summary 
before coming to the end of the program. 

In box 4, we might have found that the action item was a call 
termination (box 13). Now, we must obviously disconnect the two 
partners. At the same time, we shall also choose new action times 
for both partners (box 14). These new action times are then entered 
into the time status record (box 15). 

If we were unable to assign the call to the primary route (box 16), 
we shall update the pertinent statistic. Then, we shall try the avail- 
ability of a secondary routing (boxes 17, 18) and return to the main 
leg of the program. If the secondary routing is also unavailable for 
placement of this call (box 10), we shall update our statistics again. 
Then we shall assign to the originating subscriber randomly a new 
call initiation time (box 20). This new call initiation time is entered 
into the time status record and we return to the main leg of the 
program. 

This randomly chosen partner might turn out to be busy with 
another call (box 22). We will consider this another lost call and 
add to our statistics. Naturally, the initiator will cancel his request 
for this call. But, we must give him another chance. Therefore, we 
assign to him another random action time for call initiation (box 20). 
After entering this time in the TSR (box 21) we return to the main 
program. 

Table 1 shows the major elements which might make up the state 
vector in the time status record. The controlling element is the 
action time, ¢. Other entries include the status of 24 subscribers, 
8 switchboards, and 4 trunk lines. The subscribers are numbered 
from 1 through 24. The initiating partner is symbolically designat- 
ed I; the terminal partner T. Busyness of the switchboards and 
trunks is indicated by counting the busy channels. Therefore, 
C—o would indicate that no calls are going through switchboard C ; 
or, 23—3 would indicate that the trunk 23 is fully occupied carrying 
three conversations. 

The time status record could be completed by seven statistics to 
describe the total number of calls initiated in the system, the total 
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number of completed calls, the number of calls lost in the switch- 
boards (due to overload), the number of calls lost in switching 
centers, the number of calls going over the primary routes, and the 
number of calls going over the secondary route. 


TABLE 1. —- State vector in TSR 


Element Number Example Comments 
1. Time I t Time of last event incorporated in 
Sie 
2. Subscribers 24 1 —ÏI I, T signify that these subscribers 
2 — are busy since they are initial or 
6—T terminal parties of a circuit; 
DT indicates attemptat initiation, but 
19—T other party busy. 
3. Switchboards 8 A — 2 0, I, 2 signify the number of calls 
B —I going through board, equivalent 
C —o to local link use. 
H —2 
4. Trunks 4 I2—o 0, I, 2, 3 signify the number of cir- 
23 — 3 cuits occupied in any trunk. 
3455) 
41 — 2 
5. Statistics 7 S —I Total No. of initiated calls. 


S —2 Total No. of completed calls. 
S — 3 No. of calls lost in switchboard. 


S  —4 No. of calls lost in switching cen- 
ters. 

S —5 No. of calls lost by other party 
busy. 


S —6 No. of calls getting first route. 
No. of calls getting second route. 


| 


Total 44 


RANDOM NUMBER TECHNIQUES 


The system flow chart, Fig. 5, shows several boxes where we 
need random numbers during the simulation. For example, in 
box 5 we must randomly choose a partner for the subscriber who 
initiated the call. In this simplified model the propensity to inter- 
act is assumed the same, for any two subscribers, regardless of the 
distance which separates them. Therefore, we may make a selection 
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from uniformly distributed numbers ranging from or to 24. Usually, 
tables of random numbers are not arranged in that order and we 
shall briefly explain how to use these tables. We may take successive 
pairs of decimal random numbers (for example, from Tippett’s 
tables) and reduce them modulo 24 if they are in the range from 
00 to 95, rejecting numbers 96 through 99. Finally, we add one to 
the reduced number for the partner assignment. In mathematical 
notation, 


Pe= X,)(mi0d?24) +1 


where P, = assigned random partner, and X, = random number 
from the set 00 (OI) 05. 


For example, reading Tippett’s tables downward, the number 
sequence 23 05 I4 38 97, II ... assigns the random partners 
24 6 15 15 12 ... We reject 97 as being outside the sample range. 
Partner 15 appears twice ; we might interpret this as a call and 
simultaneous cancellation, as would be obtained when we dial a 
wrong number. 

Rectangular random numbers between 00 and 99 are also used 
to set up action or most imminent event times. We build up action 
times by first differences. For example, time increments of 64 Io 
51 60 37... in Tippett’s tables would predetermine event or 
action times of 64, 74,... Table 2 shows what action times the 
subscribers S;(¢ = I, 2, 3,..., 24) had in the single simulation run 
shown here. For example, subscriber 2 went into action at times 
65, 70, 100, 173, ... The table was terminated at the right when- 
ever a subscriber’s action time exceeded the predetermined limit of 
300 time units. Thus, the three subscribers 16, 20, and 21 went into 
action only four times, while subscriber 23 went into action eight 
times during the same interval. 

The distinction between call initiation and call termination 
(box 4) does not require the use of random numbers ; it is condi- 
tional upon the state of the system. If a subscriber is not busy and 
is called into action via the most imminent event in the TSR, the 
action becomes automatically a call initiation or an attempted call 
initiation. Conversely, if a subscriber is busy and the TSR calls for 
action, it will automatically terminate his call. Furthermore, if a 
subscriber terminates a call, his partner is automatically placed into 
the same status. A new action time is then required for both part- 
ners. These action times are always taken from Table 2, which may 
be set up in advance. 
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TABLE 2. — Action times in TSR 


S; Time 

I 89 93 IOI 232 260 325 

2 65 70 TOO 173 180 266 349 
3 8 49 135 168 229 Ce 301 
4 38 69 160 198 220 274 345 
5 73 94 149 215 245 304 

6 51 IIO 183 242 269 316 

37 18 97 184 218 273 366 

8 50 132 175 238 248 307 

9 63 161 198 236 296 322 

10 31 17 146 169 203 239 291 343 
eit 32 69 84 181 185 2 364 
12 58 IOI 17 180 278 358 

13 44 68 107 133 228 320 

14 22 109 194 219 237 308 

15 59 96 179 195 254 345 

16 55 139 176 265 327 

17 18 84 85 174 201 279" ~359 
18 18 WP 164 165 264 290 339. 
19 8 97 150 213 236 259 314 
20 ar 98 174 239 302 

an 50 10g 186 276 349 

22 49 94 141 197 - 204 348 
23 44 100 181 185 192 Big 221 235 302 
24 13 38 116 207 202 314 


In a more sophisticated communications network simulation, 
random events would be more realistically approximated by expo- 
nential holding times, Poisson distributed intervals, Erlang distri- 
butions, etc. 


SIMULATION 


For the purpose of this illustration, only one run was performed ; 
it started with an empty system. The results of this run are shown 
in Table 3, from the time of the first event (¢ = 8) to the time of 
run termination ({ == 300). 

The controlling feature in this table is the action time f in the 
first column of the table. The next two columns hold the initiating 
subscriber, S,, and his cell partner, S,. If S, (the partner being 
called) is busy, we record the fact that he does not enter the inter- 
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action by placing his call number in parentheses. The fourth co- 
lumn indicates changes in system ‘* busyness ’’. Whenever a new 
call enters the system, this column will contain a plus sign ; a call 
termination is indicated by a minus sign. In the case of an attempt- 
ed, but lost call, the plus sign will be placed in parentheses ; it 
does not really count. 

The fifth column shows the assigned routing for all completed 
calls. The next column Af indicates the advance in the time status 
record in terms of the action times shown in the first column. 

The seventh column shows the activity of the 24 system subscri- 
bers. Each call is initiated by a definite party and will be denoted 
by I; his partner, the terminal subscriber, is marked T. 
A pair of arrows extend from I and T downward to indicate the 
length of this call from initiation to termination. In case a subscri- 
ber attempts to initiate a call but cannot complete it for sundry 
reasons, the I will be crossed f. 

The next column has eight entries to record the switchboard 
activity from A to H ; 0, I, or 2 indicate the number of calls routed 
through the board. The following column records activity on the 
four trunks : 12, 23, 34, and 41 ; counting of channels from 0 to 3 
is used to designate trunk busyness. 

For the sake of expediency, not all the statistics mentioned earlier 
are shown in this table. Two columns reflect the statistics of inter- 
est. The first, S,, tells us how many calls have been placed in the 
system since start of the simulation. The second, AS,, records the 
number of live calls which go through the switchboards at any 
given time. 

Starting with an empty system, the first action time (¢ = 8) will 
obviously cause a call initiation. It involves partners 2 and 3. From 
this time on, additional calls will enter the system. Sooner or later, 
some of them will also terminate, the maximum holding time being 
99 units. Therefore, the system can be expected to settle down to a 
steady state operation. 

For the purpose of discussion, we select a few lines of this table. 
For example, the first line shows that initial action in the system 
took place at time 8. At this time subscriber 3 tried successfully 
to place a call to subscriber 2; one call entered the system. The 
required routing is simple ; the call goes only over switchboard A. 
This call terminated at time 49 by action of subscriber 3; the 
arrows in Table 3 extend to that position. This being the first call 
in the system, the switchboard activity sub-vector reads 100000 
0 o. Since no trunks are involved, the trunk sub-vector reads 000 0. 
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The pertinent statistics sub-vector is I, I. 

After some other activity, subscriber 4 successfully places a call 
to subscriber 14 at time 38. The assigned routing is B 2 3 E. There- 
fore, the switchboard sub-vector changes from 10100210 to 
LHTOT 270: the trunk. sub-vector, from 1:00 1 to 11 0 1. The 
two statistics are updated from 5 3 to6 4. At the same time, how- 
ever, subscriber 24 tried to cali subscriber 14. He gets a busysignal 
and the table entry reads f under 24. 

At time 49 the TSR indicates another action by subscriber 3 
(Table 2). Since he is busily engaged in conversation with subscri- 
ber 2, the action will mean termination of this call. Thus, the fourth 
column in Table 3 contains a minus sign. The board sub-vector 
decreases in activity from 11102210to0r1102210. Trunk 
lines are not involved and the trunk sub-vector remains unchanged. 

In summary, Table 3 portrays the activity which permeates the 
system during the simulation run. Naturally, we must summarize 
this information in some way so as to obtain a better over-all view 
of the system and its operation. The dynamic system response is 
shown graphically in Fig. 6. The horizontal scale shows time elapsed 
in basic units (say, minutes). The vertical scale indicates counts or 
frequencies. 

The first curve shown is S,. This step function portrays cumula- 
tive system activity. Because of the simplified assumptions made, 
the curve has a linear upward trend. The second curve AS, shows 
the number of calls which are in existence in the system at any one 
given time, the system ‘‘ busyness ’’. Naturally, this too is a fluc- 
tuating step function. Ignoring the initial transient, it attains soon 
a steady state average of approximately 5 (calls). During the first 
300 minutes of simulated time, the system’s busyness drops once to 
two calls (at time 184) ; it also attains a maximum of 7 simultaneous 
calls at times 116, 160, 237, 265, and 269. We observe in passing 
that the system is capable of carrying eight simultaneous calls, at 
most. 

Fig. 6 indicates that steady state commences approximately at 
time 50. As we have pointed out earlier, one cannot determine 
‘ precisely ’’ a time of transition from transient to steady state, 
even in a system as simple as this one. Nevertheless, we can state 
that steady state has started after time 50, with high probability. 
On the other hand, the probability that the system attains steady 
state before time 40 is considerably smaller. Finally, trend ana- 
lysis, or some other statistical tool, would show that the system 
“ certainly ” is not in steady state between times 0 and 30. 
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EXTENSION PROBLEMS 


Thus far we have considered the problem of simulation in general, 
and in particular, the simulation of communications networks. The 
simple demonstration of the preceding section should suffice to 
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illustrate the power of digital simulation programs. However, before 
closing this expose, we might just try to take a more searching 
look at the general problem, with a view toward its limitations. 

The first point we want to bring up is the distinction between 
analysis and synthesis. Generally speaking, simulation is a tool of 
analysis. Given a real system, or only its concept, it can be simu- 
lated to any desired degree of truthfulness on digital computers. 
We need not concern ourselves with the fact that some simulations 
may use up excessive amounts of computer time. If such simulations 
become necessary, man will devise faster computers (1). Therefore, 
as an analytical technique, simulation can provide a maximum of 
information about a real, hypothesized, or postulated system. In 
fact, the analyst often does not even know what to look for — 
but surely finds IT during computer runs. 

This is then the crux of the matter. Simulation is an excellent 
analytical tool ; but it provides little, if any, information about the 
inverse problem of synthesis. 

Consider a simple, but typical example. The computer analyst 
is asked to synthesize the ‘* best ’’ subroutine for the computation 
of y = sin x. The word “ best ’’ presents an immediate stumbling 
block. What is meant by * best ”’ ? The fastest routine, regardless 
of memory space used ; or a routine requiring least memory space ; 
or some sort of a weighted compromise ? In the latter case, who 
provides the weighting factors for the compromise ? 

Now suppose that we want the fastest routine, regardless of 
memory space requirements ; also, that the computer and its 
vocabulary are specified. Present-day custom is to assign this job 
to a programmer. He sits down and does some thinking, like hund- 
reds before him, and comes up with a flow chart, and a piece of 
coding. No sooner than he turns his results over to the ‘* boss ”’ 
that somebody takes a look at them and finds a better way of 
doing the whole job. This is known in programmer’s parlance as 
‘ the second round ’’. And so it goes, somebody improves this 
program, round after round. Occasionally, the improvement is a 
real “ breakthrough ”’ by virtue of a radically novel scheme. How- 
ever, it, too, will eventually make room for a better one. Thus, no 
programmer can ‘* guarantee ”’ that his is the final, ultimate defini- 
tive, THE subroutine. What we do, is to test each subroutine for 


(1) Such was indeed the case in the design of LARC (by Remington Rand UNI- 
VAC) for the Atomic Energy Commission whose simulation requirements exceeded 
existing computer capacities (1957) by a factor of one-hundred. 
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speed and to record the results and any improvements. But the 
problem of synthesis remains unsolved ; only exhaustive analysis 
has been performed on each routine. As a result, all but one have 
yielded to their respectively superior successors by going unglo- 
riously to the circular file. 

We see now that simulation leads to extensive testing. We inter- 
pret this analytical phase as synthesis by trial and error. But a real 
synthesis has not been accomplished. In fact, the reason that so 
much work is done on computers is often traced to their inability to 
synthesize ; or to man’s limited ways in the search for a rationale 
leading to complete syntheses. ‘* Comparative analysis ? would be 
a better term in almost all cases where ‘‘ optimum ”’ results by 
synthesis are being claimed. 

Nevertheless, we are in a position to propose a rationale for 
synthesis. However, in the same breath we shall show its nonfeasi- 
bility. Computers have a finite vocabulary. Therefore, we could 
write a random program generator, e. g., a program which produces 
random sequences of instructions. These sequences could be tested 
to determine their ‘ output ’’. One of the myriad routines would 
be the desired ‘* best ’’ subroutine for y = sin x. This is the elec- 
tronic analog to the old story of the monkey hitting random keys on 
a typewriter ; ultimately he will produce not only the works of 
Shakespeare, but also all the as yet unwritten great epics of the 
future. Random instruction sequence testing, too, might take 
Googolplex years. Thus we must reject this rationale on realistic 
grounds and revert to limited synthesis by comparative analysis. 

An extension of the synthesis problem is the augmentation 
problem, e.g., improvement by synthesis. To illustrate, let us 
return to the simulation of communications networks. Suppose a 
very good simulation program has been written for an existing or an 
hypothesized network. Now we test the system on the computer 
by feeding it inputs with realistic traffic rates, holding times, etc. 
Suppose we find that the system response indicates saturation in one 
or more small areas or subsystems. For example, think of a global 
communications network which might reveal under test that the 
Canadian area could become clogged, or saturated, while the rest 
of the system is operating normally. The design engineer would 
now like to ‘‘ augment ’’ the system. He wants to remove the 
danger of saturation by adding more stations, or additional equip- 
ment in the saturated area(s) only. In the case of additional stations, 
the question is simply where to place them to produce the best 
results. This is a simple question ; but we are hard put to find an 
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answer. Simulation techniques may be used to provide synthesis 
by repetitive analysis. Again, this is a trial-and-error method and 
true synthesis is now not possible, although the method of exhaus- 
tion would provide a theoretical, but unfeasible solution. Limited 
exhaustion through a combination of Monte Carlo and analytical 
techniques will probably remain in use for some time to come. 

Let us finally look at routing doctrines. There exist several con- 
ceptually different techniques and as many variations. Now, if we 
are given a communication system design AND one stated doctrine, 
it can be tested and analyzed by simulation. A second doctrine or 
several variants may also be tested and compared with the first by 
subjecting them to computer experimentation. Finally, one doctrine 
will be selected as having the greatest merits, or optimum utility 
in some sense. Yet it is clear that comparative testing of additional, 
new routing doctrines could continue almost forever. However, this 
work might never produce any spectacular results in the sense that 
an ‘optimum ’’ routing doctrine is discovered. 

Relevant Monte Carlo techniques for iterative augmentation are 
also of interest. They may be written so as to show on the monitor 
printer results only if the most recent run has bested all preceding 
ones, especially the last best one (which in turn had bested all its 
predecessors). The observer of such a program gets a really good 
feeling for the incremental gains attained. Initial printouts follow 
one another rapidly. Soon, the flow slows down to a point where 
the intervals between printouts grow longer and longer. Eventually, 
a decision must be made whether or not to continue the program ; 
this decision may also be programmed. But, as in the case of the 
ruined gambler, one is left with that empty feeling that given just 
one more chance, the results might have been of dramatic import. 

Many extension problems which give rise to simulation programs 
can be traced to a need for synthesis by analysis. Therefore, we 
may surmise that a good deal of the work done on digital compu- 
ters is done for that reason by iterative analysis testing. Small 
wonder then that the computer industry flourishes ! It, too, must 
solve many of its own synthesis problems to satisfy demands for 
new hardware and more versatile software.’ 


SUMMARY AND CONCLUSIONS 


This has been an attempt at providing an expose of a specific 
problem area which has attracted considerable attention. We have 
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delved into some of the complex issues involved in digital computer 
simulation ; we have even given a simple example to illustrate 
applicable techniques. Let us conclude by admitting that any 
expose of such a vast topic must remain necessarily superficial. 
In fact, as we survey the field, we observe that even the experts 
continually discover a wealth of new material. Attempts at nation- 
wide coordination of these efforts have been made by many societies 
and organizations. Most notable in this connection are the meetings 
held at Ann Arbor [10, 11] and the efforts made by the Operations 
Research Society. 

Computers were originally introduced as labor saving devices. 
Now, ten years later, new industries have sprung up, new job cate- 
gories and new problems are facing us. Where the original idea had 
been to cut down the workweek our thinking had to be revised. 
While there will be more time available to large sectors of the labor 
force, this time must be invested in education and self-improve- 
ment. Only this process will give everybody access to new and 
better jobs. More important, it will broaden our horizons and provide 
greater challenges. In that sense, let us hope that the electronic 
revolution which began ten years ago with computers may ulti- 
mately lead towards the highest educational standards for our 
whole society. 


REFERENCES 


[1] Arkin, H. and Cotton, KR. R., Tables for statisticians. Barnes & Noble, 
1950. 

2] CHURCHMAN, C. W., ACKOFF, R. L., ARNOFF, E. L., Introduction to 
operations vesearch. John Wiley, New York, 1950. 

[3] GRABBE, E. M., Ramo, S. and WooLDRiDGE, D. E., Handbook of auto- 
mation, computation, and control (Volumes 1 & 2). John Wiley, New 
York, 1958. 

[4] GREGORY, R. H. and Van Horn, R. L., Automatic data-processing 
systems. Wadsworth Publishing Company, 1961. 

[5] Hastines, C. Jr., Approximations for digital computers. Princeton Uni- 
versity Press, 1955. 

[6] IkLÉ, F. C. and Hammer, C., Intercity telephone and airline traffic 
velated to distance and the propensity to interact. Sociometry, December 
1957. 

[7] LeGAULT, R. R., MILLER, J. M. and GoLDBERG, M., Four papers con- 
cerning the modeling of communications systems. University of Michi- 
gan Report 4592-1-T, June 1961. 


[8] 
[9] 
[10] 


[11] 
[12] 


[13] 


[14] 


COMPUTERS AND SIMULATION 249 


McKinsey, J. C. C., Introduction to the theory of games. McGraw-Hill, 
New York, 1952. 

NEUMANN (von) J. and MORGENSTERN, O., Theory of games and economic 
behavior. Princeton University Press, 1953. 

Proceedings, Symposium on digital simulation techniques for predicting 
the performance of large-scale systems, 23-25 May 1960. Ann Arbor, 
Michigan (Report 2354-33-X). 

Proceedings, Thivd war games symposium, 6-8 October 1960. Ann Arbor, 
Michigan (Report 36943-18-X). 

RALSTON, A. and Wir, H. S., Mathematical methods for digital com- 
puters. John Wiley, New York, 1960. 

REITWIESNER, G. W., An ENIAC determination of 7 and e to more than 
2000 decimal places. Mathematical tables and other aids to compu- 
tation, January 1950. 

UNESCO, Proceedings of the International Conference on Information 
Processing, Paris, 15-20 June 1959. R. Oldenbourg (München), 1959. 


[15] Vazsony1, A., Scientific programming in business and industry. John 


Wiley, New York 1958. 


CYBERNETICA (Namur) 
Vol. IV — N° 4— 1961 


